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Abstract 

The production of dijets in diffractive deep inelastic scattering has been measured 
with the ZEUS detector at HERA using an integrated luminosity of 61 pb _1 . The 
dijet cross section has been measured for virtualities of the exchanged virtual 
photon, 5 < Q 2 < 100 GeV 2 , and 7*p centre-of-mass energies, 100 < W < 250 
GeV. The jets, identified using the inclusive kr algorithm in the 7*p frame, 
were required to have a transverse energy E^^ et > 4 GeV and the jet with the 
highest transverse energy was required to have E^- Ct > 5 GeV. All jets were 
required to be in the pseudorapidity range —3.5 < ?7j* ct < 0. The differential 
cross sections are compared to leading-order predictions and next-to-leading- 
order QCD calculations based on recent diffractive parton densities extracted 
from inclusive diffractive deep inelastic scattering data. 
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1 Introduction 



Diffractive events in deep inelastic scattering (DIS) are characterised by the presence of a 
fast forward proton, a large rapidity gap (LRG) - an angular region between the scattered 
proton and the dissociated photon with no particle flow [1-6] - and a dissociated virtual 
photon 7*. In recent years perturbative QCD (pQCD) has become a successful tool for 
describing diffractive events [4-7]. The cross section for diffractive DIS processes can be 
described by a convolution of universal diffractive parton distribution functions (dPDFs) 
and process-dependent coefficients, which can be calculated in pQCD [8]. At HERA, 
dPDFs have been determined using inclusive diffractive DIS data [4-6] . 

This paper presents measurements of dijet production in diffractive neutral current DIS 
with the ZEUS detector at HERA. The presence of a hard scale in such a process, either 
the virtuality of the photon or the large jet transverse momentum, is well suited for 
a pQCD analysis. Dijet processes are particularly sensitive to the density of gluons in 
the diffractive exchange (i.e. via 7* g — > qq, as shown in Fig. [1]), and gluons have been 
shown to carry most of the momentum of the colourless exchange [4,5,9]. The measured 
differential cross sections are compared with leading-order (LO) and next-to-leading-order 
(NLO) QCD predictions using the available dPDFs. The results presented here benefit 
from higher statistics compared to previous measurements of the same process [10]. 



2 Experimental set-up 

This analysis is based on 61 pb -1 of data collected with the ZEUS detector at the HERA 
collider during the 1999-2000 data-taking period. During this period, HERA collided 
either electrons or positrons^ of 27.5 GeV with protons of 920 GeV at a centre-of-mass 
energy of -y/i = 318 GeV. 

A detailed description of the ZEUS detector can be found elsewhere [11]. A brief outline 
of the components that are most relevant for this analysis is given below. 

Charged particles are tracked in the central tracking detector (CTD) [12], which operates 
in a magnetic field of 1.43 T provided by a thin superconducting coil. The CTD consists 
of 72 cylindrical drift chamber layers, organised in 9 superlayers covering the polar-angle 
regional 15° < 9 < 164°. The transverse-momentum resolution for full-length tracks is 

1 In the following, for simplicity, the word positron will be used to denote both electrons and positrons. 
The integrated luminosity for e~p data is 3 pb _1 , while for e + p data is 58 pb 1 . 

2 The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. 
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o{Vt)Ivt = 0.0058p T © 0.0065 © 0.0014/p T , with p T in GeV. 

The high-resolution uranium-scintillator calorimeter (CAL) [13] consists of three parts: 
the forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. Each part 
is subdivided transversely into towers and longitudinally into one electromagnetic sec- 
tion and either one (in RCAL) or two (in BCAL and FCAL) hadronic sections. The 
smallest subdivision of the calorimeter is called a cell. The CAL energy resolutions, 
as measured under test-beam conditions, are a(E)/E = 0.18/V^E for electrons and 
a(E)/E = 0.35/ VE for hadrons, with E in GeV. 

During the 1999-2000 data-taking period, the forward plug calorimeter (FPC) [14], located 
in the beam hole of FCAL, extended the pseudorapidity coverage of the calorimeter up 
to values of rj ~ 5. It consisted of a lead-scintillator calorimeter read out by wavelength 
shifters and photomultipliers. 

In order to improve the detection of positrons scattered at low angles, the angular coverage 
in the rear direction was extended by means of the small rear tracking detector (SRTD) 
[15, 16]. The SRTD consists of two planes of 1 cm wide and 0.5 cm thick scintillator 
strips glued on the front of RCAL. The orientations of the strips in the two planes are 
orthogonal. Scattered positrons were also detected in the rear hadron-electron separator 
(RHES) [17], a matrix of more than 10000 silicon diodes 400 yim thick inserted in the 
RCAL. 

The luminosity was measured using the bremsstrahlung process ep — > epy with the 
luminosity monitor [18], a lead-scintillator calorimeter placed in the HERA tunnel at 
Z = -107 m. 

3 Kinematics 

Dijet production in diffractive DIS (ep — > e + p + jl+j2 + X') is characterised by the 
simultaneous presence of a scattered positron, a scattered proton p that escapes unde- 
tected down the beam pipe, and the photon-dissociative system X, which contains the 
dijet system jl + j2, produced in the hard scattering along with the rest of the hadronic 
system X' (see Fig. [T]). Deep inelastic scattering of a positron on a proton is described by 
the following kinematic variables: 

• s = (P + k) 2 , the squared ep centre-of-mass energy, where P and k indicate the 
incoming proton and the incoming positron four-momenta, respectively; 

• Q 2 = —q 2 = —{k — k') 2 , the virtuality of 7*, where kl is the four-momentum of the 
scattered positron; 

• W 2 = (P + q) 2 , the centre-of-mass energy squared of the 7*p system. 
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Diffractive events are further characterised by the variables: 

• Mx, the invariant mass of the photon-dissociative system; 

• t = (P — P') 2 , the squared four- momentum transfer at the proton vertex, where P' 
denotes the four-momentum of the scattered proton; 

• x P = (P — P') ■ q/P ■ q, the momentum fraction lost by the proton; 

• (3 = Q 2 /2(P — P') ■ q, a measure of the fractional momentum of the diffractive exchange 
carried by the struck parton. 

The description of the dijet system in the hadronic final state requires the use of additional 
variables: 

• Zp, the fraction of the momentum of the diffractive exchange carried by the parton 
participating in the hard process and defined as 

_ q-v 

where v is the four-momentum of the parton originating from the diffractive exchange; 

• sc 7 , the fractional momentum of the virtual photon participating in the hard process. 
In DIS, x 7 is expected to be unity (direct photon). However, some models introduce 
the concept of a resolved virtual photon, where the 7* can fluctuate into a partonic 
state before participating in the hard interaction. For resolved photon processes, x 1 
is expected to be lower than unity. The variable x 1 is defined as 

P-u 



P-q' 

where u is the four-momentum of the parton originating from the virtual photon. 



(2) 



4 Theoretical models 



4.1 QCD factorisation in diffraction 

The cross section for diffractive DIS processes at fixed s depends in general on four 
independent variables, which are usually chosen to be Q 2 ,j3,x F and t. According to the 
QCD factorisation theorem [8], the cross section for inclusive diffraction, a(j*P ~~ * Xp), 
can be written as 
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This expression is valid at fixed x F and t and for scales sufficiently large to permit the 
use of pQCD. The sum runs over all partons. The partonic cross-section cx 7 **(Q 2 , 0) for 
the hard subprocess involving the virtual photon and the parton i is calculable in pQCD. 
The functions fP(x P , t, /3, Q 2 ) are the dPDFs: they describe the probability to find in 
the proton a parton of kind i carrying a fraction x F ■ [3 of its momentum with a probe 
of resolution Q 2 under the condition that the proton stays intact, with a momentum loss 
quantified by x P and t. For diffractive production of dijets, Eq. ([3]) is rewritten as: 

i 

^ ^ f dQ 2 [ d£ aTMQ 2 , ^)f t D (xjp, t, f , Q 2 ), 



dx~dt~ ^ I u " / ^ 
t=i,q,9 



where now z P is the variable sensitive to the dPDFs and the subprocess cross section a 1 " 1 
is replaced by the cross section, cj. 1 , for the reaction 7*2 — > jl j2. 

At HERA, the dPDFs have been determined within the QCD DGLAP formalism [19-22] 
by means of fits to inclusive diffractive DIS measurements with a procedure similar to that 
used to extract the standard proton PDFs from inclusive DIS data [23-28]. Consistency 
between the measured cross sections for semi-inclusive processes and calculations using 
these dPDFs represents an experimental proof of the validity of the QCD factorization 
hypothesis in diffraction [10,29]. 

Most of the dPDF parameterisations use Regge phenomenology arguments [30] to factorise 
the (x P ,t) from the (fi,Q 2 ) dependence. In the Regge approach, diffractive scattering 
proceeds via the exchange of the Pomeron trajectory. The dPDFs are then written as 
the product of the Pomeron flux (dependent on Xp and t) and parton distributions in 
the Pomeron (dependent on (3 and Q 2 ). For x P values substantially larger than 0.01, the 
contribution of the subleading Reggeon trajectories may also have to be added. 



4.2 NLO calculation 

Predictions for diffractive dijet differential cross sections were calculated at order a 2 s 
with the program Disent [31] adapted for diffractive processes. The calculations were 
performed in the MS renormalisation scheme with five active flavours and the value of 
the strong coupling constant set to a s (M z ) = 0.118. The predictions were obtained 
with the renormalisation scale, /ir, equal to E^^, where -Exji i s the transverse energy 
of the highest transverse energy jet in the event (the leading jet) as measured in the 7*p 
centre-of-mass frame. The factorisation scale was set to Q 2 . 

The following dPDFs were used: 
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• the ZEUS LPS+charm [6] - the result of an NLO DGLAP QCD fit to the inclusive 
diffractive structure functions measured by the ZEUS experiment with the leading 
proton spectrometer (LPS). In order to better constrain the dPDFs, measurements 
of D* production cross section in diffractive DIS [32] were also included. The fit was 
restricted to the region x P < 0.01; 

• the HI 2006 dPDFs [5] - the result of an NLO DGLAP QCD fit to a sample of inclusive 
diffractive structure functions measured by the HI Collaboration. Two different pa- 
rameterisations are available (Fit A and B) which differ in the gluon distribution. The 
fit was restricted to the region Q 2 > 8.5 GeV 2 , z P < 0.8. Since the HI measurements 
were not corrected for the contribution due to events where the proton dissociated 
into a low-mass state, in the comparison the calculations were renormalised by a fac- 
tor 0.87 [5]; 

• the Martin- Ryskin- Watt 2006 (MRW 2006) dPDFs [7] - the result of a fit to the same 
data set as for the HI 2006 fit. Regge factorisation is assumed only at the input scale. 
The dPDFs are then evolved with an inhomogeneous evolution equation analogous 
to that for the photon PDFs. The inhomogenous term accounts for the perturbative 
Pomeron-to-parton splitting. 

The only theoretical source of uncertainty considered was that coming from the NLO 
calculations. This uncertainty was estimated by varying /xr by factors of 0.5 and 2. 
Uncertainties of more than 20% were obtained. To compare with the data, the NLO 
predictions at the parton level were corrected to the hadron level using factors extracted 
from a MC program (see Section [5]). The corrections were typically of the order of 10%. 

5 Monte Carlo simulation 

Monte Carlo (MC) simulations were used to correct the data for acceptance and detector 
effects. Two different MC generators were used, Rapgap [33] and Satrap [34]. 

The Rapgap MC is based on the factorised-Pomeron approach. The events were gen- 
erated using the HI fit 2 dPDFs [4]. No Reggeon contribution was included in this 
simulation. The parton-shower simulation is based on the Meps [35] model. Resolved 
photon processes were also generated using Rapgap with the GRV-G-HO [36] photon 
PDFs. Since the relative contributions of direct and resolved photon processes to the 
total cross section are a priori unknown, the Rapgap direct and resolved samples were 
weighted in order to best describe the data. The Rapgap MC was also used to extract 
the hadronisation corrections for the NLO calculation. 

Satrap is based on the Golec-Biernat-Wiisthoff (GBW) dipole model [34] and is inter- 
faced to the Rapgap framework. The parton-shower simulation in Satrap is based on 
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the Colour Dipole Model (CDM) [37]. This MC does not include the resolved-photon 
contribution to the 7*p cross section. 

To estimate the inclusive DIS background, a sample of events was generated with Djan- 
GOH [38]. 

All the above MC programs are interfaced to the Heracles [39] event generator for the 
simulation of QED radiative processes and to Jetset [40] for the simulation of hadronisa- 
tion according to the Lund model [41]. QED radiative corrections were typically between 

5 and 10%. 

The ZEUS detector response was simulated with a program based on Geant 3.13 [42]. 
The generated events were passed through the detector simulation, subjected to the same 
trigger requirements as the data, and processed by the same reconstruction and offline 
programs. The average of the acceptance-correction values obtained with Rapgap and 
Satrap was used to correct the data to the hadron level. 

6 Event reconstruction and data selection 
6.1 DIS selection 

A three-level trigger system was used to select events online [11,43]. In the third- level 
trigger, a DIS positron candidate and energy deposition in the FPC lower than 20 GeV 
were required. The scattered positron was identified both online and offline using a neural- 
network algorithm [44] . The reconstruction of the scattered positron variables was carried 
out by combining the information from CAL, SRTD and HES. In order to select a DIS 
sample the following requirements were applied [45,46]: 

• the positron found in the RCAL had to lie outside a rectangular area of size [-14, +12] 
cm in X and [-12, +12] cm in Y, centred around the beam pipe. Further cuts on the 
fiducial area of the impact point of the positron on the RCAL surface were applied in 
order to exclude regions with significant inactive material [47]; 

• the energy of the scattered positron had to be greater than 10 GeV; 

• the vertex of the event had to be in the range |Zvtx| < 50 cm to reject non-ep back- 
ground. 

The four-momentum of the hadronic final-state X was reconstructed using energy-flow 
objects (EFOs), which combine the information from the CAL and the CTD [48]. The 
EFOs were corrected for energy losses due to the inactive material present in the detector 
[49]. The variable 5 = J2i= e efo(^ ~Pz,i), where the sum runs over the scattered positron 
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and all the EFOs, was required to be 45 < 5 < 65 GeV. The variables E( and pz,i denote 
the energy and the Z-component of the momentum of each term of the sum. 

The Q 2 and W variables were determined using the double-angle method [50]. Events 
were accepted if 5 < Q 2 < 100 GeV 2 and 100 < W < 250 GeV. 

6.2 Jet selection 

The fc^-cluster algorithm in the longitudinal invariant mode [51] was applied to the cor- 
rected EFOs in the photon-proton centre-of-mass system (^*p frame) to reconstruct the 
jets. The jet variables in the 7*p frame are denoted by a star. After reconstructing the 
jets, the massless four-momenta were boosted to the laboratory frame where further en- 
ergy corrections were determined and propagated back into the transverse energy of the 
jet, Ejijet- Such corrections, obtained from a MC study, improved the correlation between 
hadron- and detector- level transverse energy of the jets [45]. The dijet sample was defined 
by requiring the events with at least two jets to fulfill the following constraints: 

• > 5 GeV and E^ - 2 > 4 GeV, where the labels jl and j2 refer to the jets with the 
highest and the second highest transverse energy, respectively; 

• —3.5 < r]* ct < 0, where rj* et is the pseudorapidity of any of the jets; 

• the pseudorapidity of the selected jets, boosted to the laboratory frame, had to lie in 
the range |f7j^ B | < 2. 

6.3 Diffract ive selection 

Diffractive events are characterised by low values of x P and by the presence of a LRG. 
The following selection criteria were applied [45,46]: 

• E-ppc < 1 GeV, where -Efpc is the total energy in the FPC. The requirement of activity 
compatible with the noise level in the angular region covered by the FPC is equivalent 
to a rapidity-gap selection; 

• x°p S < 0.03 where x^ s is the reconstructed value of x P and is defined as: 

r obs = Q 2 + Mj 
Q 2 + W 2 ' 

The mass of the diffractive system, Mx, was reconstructed from the EFOs. The 
cut on x°p S reduces the contribution of Reggeon exchange and other non-diffractive 
background. 
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After these cuts, the selected sample is still contaminated by diffractive events in which the 
p dissociated into a low-mass system. This contamination was estimated by MC studies 
to be /pdiss = (16 ± 4)% [52] and was subtracted from the measurements independent of 
the kinematics. 

The contamination of the non-diffractive background as a function of the applied diffrac- 
tive selection cuts is shown in Fig. [2j through the distribution of t?max, where 7/max is 
the pseudorapidity in the laboratory frame of the most forward EFO with energy higher 
than 400 MeV, before and after applying cuts on the Eppc and on x^ s . The disagreement 
between the measured and the simulated distributions is the reason for not applying any 
requirement on 7/max, as was done in previous analyses [9,32,53]. After the -Efpc and 
x^ s cuts, the non-diffractive background from Djangoh was estimated to be 2.4% of 
the total selected events and neglected in further analysis. After all cuts, 5539 events 
remained. 

7 Systematic uncertainties 

The systematic uncertainties of the measured cross sections were calculated by varying 
the cuts and the analysis procedure. The systematic checks were the following: 

• the energy measured by the CAL was varied by ±3% in the MC to take into account the 
uncertainty on the CAL calibration, giving one of the largest uncertainties. Deviations 
from nominal cross section values were of the order of ±5%, but reached ~ 15% in 
some bins; 

• the energy scale of the scattered positron was varied in the MC by its uncertainty, 
±2%. The resulting variation of the cross sections was always below ±3%; 

• the position of the SRTD was changed in the MC by ±2 mm in all directions to account 
for the uncertainty on its alignment. The change along the Z direction gave the largest 
effect and in a few bins caused a cross section variation of ±2%; 

• the model dependence of the acceptance corrections was estimated by using separately 
Rapgap and Satrap for unfolding the data. The variations from the central value 
(obtained using the average between Rapgap and Satrap) were typically of the order 
of ±5% but reached ~ ±10% in some bins. 

The above systematic uncertainties, except those related to the energy scale of the calorime- 
ter, were added in quadrature to determine the total systematic uncertainty. The un- 
certainties due to the energy scale and the proton dissociation subtraction (±4%) were 
added in quadrature and treated as correlated systematics. The energy scale uncertainty 
is quoted separately in the tables. 
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The stability of the measurement was checked by varying the selection cuts as follows: 

• the cut on the FPC energy was varied by ±100 MeV in the MC; 

• the cut on the scattered-positron energy was lowered from 10 to 8 GeV; 

• the fiducial region for the positron selection was enlarged and reduced by 0.5 cm; 

• the lower cut on 5 was changed from 45 to 43 GeV. 

The variations of the cross section induced by these stability checks were small, within 
±2%, and were added in quadrature to the total systematic uncertainty. The uncertainty 
on the luminosity measurement (2.25%) was not included. 

The measurement was repeated with the addition of a cut on the value of //max- This 
estimates the uncertainty on the purity of the diffractive selection. A cut of ?7max < 2.8 
was applied. The cross sections increased by ~ 5% and the change was concentrated at 
high values of x^. No significant dependence on other variables was observed. This vari- 
ation is listed in the tables for completeness but not included in the quoted uncertainties 
of the measurement. 



8 Results and discussion 

The single- and double-differential cross sections for the production of dijets in diffractive 
DIS have been measured for 5 < Q 2 < 100 GeV 2 , 100 < W < 250 GeV and x P < 0.03, for 
jets in the pseudorapidity region —3.5 < r/* ct < 0, with E^^ > 5 GeV and E^- 2 > 4 GeV. 
The cross sections refer to jets of hadrons and are corrected for QED effects. 

The measured total cross section (given in Table [T|) is: 



a(ep -> ep + jl +j2 + X') = 89.7 ± 1.2(stat) +^(syst.) ^(corr.) pb. 

The values of the differential cross sections are averaged over the bin in which they are 
presented. For any variable k, the cross section was determined as 

do_ r N D (1- f pdiss ) 

dn CAk ' K) 

where Nd is the number of data events in a bin, C includes the effects of the acceptance 
and the QED correction factor as determined from MC, £ is the integrated luminosity 
and Ak is the bin width. 

The differential cross sections were measured as a function of Q 2 , W, x°p S , j3, Mx, E^j, 
r]j, z° bs and x° bs . The variable E^ 3 (r/j) stands for both E^^ (77^) and E^ - 2 (r/* 2 ) - in the 
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corresponding cross section, it thus contributes two entries per event. The variable 
is an estimator of z P and is calculated as 



^ _ Q 2 + W 



g 2 + M|' 



where My is the invariant mass of the dijet system. The estimator of x 7 , x° bs , is 



^AB e -^ AB +jE LAB e - 



_„LAB 

X 7 



hadr 

where the sum in the denominator runs over all the hadrons. The values of the differential 
cross sections are presented in Tables I21TTT1 and shown in Figs. [3] and HI 



8.1 Comparison to Monte Carlo models 

The Rapgap and Satrap MC programs are compared to the measured cross sections in 
Figs.[3]and|U Since the MC predictions are not expected to describe the normalisation, the 
cross sections predicted by both MCs were normalised to the data. The total correlated 
uncertainty is shown as a shaded band in the figures. The comparison with MC predictions 
shows in general a reasonably good agreement with the shape of the data. The E^ 3 
distribution is a steeply falling function as expected in pQCD (Fig. and the jets tend 
to populate the 7* fragmentation region. 

The most prominent features of the data are the rise of the cross section with x^, the 
peak at z° bs ~ 0.3 and the tail of the cross section at low x° bs values. The requirement 
of two jets with high suppresses the contribution of low values of x°p S . The relatively 
low value of the peak position in the z^ s distribution indicates that in the majority of the 
events the dijet system is accompanied by additional hadronic activity. A disagreement 
between data and Rapgap is observed at high ^ bs . In the high 2^ bs region, Rapgap 
underestimates the number of events while Satrap agrees with the data, possibly because 
of the presence of a mechanism for exclusive direct production. Most of the events are 
produced at large expected in DIS. At low x° bs , the description by Rapgap is 

improved by the addition of the resolved photon contribution (16%). 
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8.2 Comparison to NLO QCD predictions 



In Table HJ the four NLO predictions described in Section 14.21 are compared to the mea- 
sured total cross section. The central values of the predictions using the HI 2006 — FitB 
and MRW 2006 dPDFs give the best description, while those using the HI 2006 — FitA 
and the ZEUS LPS+charm dPDFs are higher in normalisation. 

The NLO predictions for the differential cross section are compared to the data in Figs. [5] 
and [6j The estimated theoretical uncertainties are shown only for the calculations using 
the ZEUS LPS+charm dPDFs and are similar for all the other calculations. For ease 
of comparison the ratios of data to the ZEUS LPS+charm prediction are presented in 
Figs. [7] and [SJ The variation due to the choice of the dPDFs is displayed with respect 
to the ZEUS LPS+charm in the same figure. In general the shape of the measured 
cross section is described by the NLO calculations within the theoretical uncertainties. 
However, only the predictions using the HI 2006 — FitB and MRW 2006 dPDFs are able 
to describe satisfactorily the data over the entire kinematic range. 

The NLO predictions for the differential cross section are compared to the data in Figs. M 
anddOl where the z° bs distribution is shown for different regions of -Epji and Q 2 . Within 
the theoretical uncertainties, the HI 2006 — FitB and MRW 2006 dPDFs are compatible 
with the data. Since the major difference between the HI 2006 — FitB and Fit A is in the 
gluon dPDF, these data have a significant potential to further constrain the gluon dPDF. 



9 Conclusions 



The single- and double-differential cross sections for the production of dijets in diffractive 
DIS have been measured with the ZEUS detector in the kinematic region 5 < Q 2 < 100 
GeV 2 , 100 < W < 250 GeV and x P < 0.03, requiring at least two jets with Sf, • t > 4 GeV 
in the pseudorapidity region —3.5 < ?7j* ct < 0.0 and the highest jet with > 5 GeV. 

Two leading-logarithm parton-shower models, Rapgap and Satrap, describe the shape 
of the measured cross sections well. The measured cross sections are able to discriminate 
between NLO QCD calculations based on different dPDFs, showing a satisfactory agree- 
ment with the calculations using the HI 2006 - FitB and MRW 2006 dPDFs. This lends 
further support to the validity of QCD factorisation in hard diffractive scattering. Since 
the dPDFs used differ mostly in the gluon content, these data may have a significant 
potential to constrain the diffractive gluon distribution. 
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a 


"stat 


^syst 


Ses 


^theor 


Adiffr 




(pb) 


( P b) 


( P b) 


(pb) 


( P b) 


(pb) 


Data 


89.7 


1.2 


+3.2 
-5.3 


+5.1 
-3.7 




+4.0 


ZEUS LPS+charm 


120.3 








+29.4 
-18.3 




HI 2006 - Fit A 


130.2 








+31.2 
-19.9 




HI 2006 - Fit B 


102.5 








+24.7 
-15.6 




MRW 2006 


99.3 








+23.4 
-14.7 





Table 1: Total cross section for the production of diffractive dijets compared to ex- 
pectations of NLO calculations using various dPDFs as indicated in the Table. The 
cross section is given for jets with -Eyji > 5 GeV , E^ 2 > 4 GeV, —3.5 < r]* ct < 
and in the range of 5 < Q 2 < 100 GeV 2 , 100 < W < 250 GeV and x P < 0.03. The 
statistical, 5 stat , uncorrected systematic, S syst , and energy scale uncertainties, 5es, 
are quoted separately. The theoretical uncertainty on the NLO calculations, <5theor; 
is quoted in the sixth column. The difference with the measured cross section with 
and without t/max cut, Adiffr, is presented in the last column. The uncertain- 
ties on the proton dissociation subtraction and the luminosity measurement are not 
presented in the table. 
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Q 2 bin 


da/dQ 2 


"stat 


"syst 




^DIFFR 


^ vjc V J 


(r>h/ dfAf 2 \ 
^pu/ Vjrbv J 




(r,h/ dp\f 2 \ 
IP D / <-* ev ) 


(r>h/ dp\l 2 \ 
^pL)/ vjt:V J 




5, 8 


7.4 


± 0.3 


+0.3 
-0.5 


+0.5 
-0.5 


0.1 


8, 12 


4.2 


± 0.2 


+0.2 
-0.3 


+0.3 
-0.3 


0.1 


12, 17 


2.6 


± 0.1 


+0.2 
-0.2 


+0.2 
-0.2 


0.2 


17, 25 


1.38 


± 0.06 


+0.09 
-0.13 


+0.08 
-0.08 


0.06 


25, 35 


0.94 


± 0.04 


+0.07 
-0.07 


+0.06 
-0.05 


0.06 


35, 50 


0.53 


± 0.03 


+0.02 
-0.03 


+0.03 
-0.03 


0.01 


50, 70 


0.27 


± 0.02 


+0.02 
-0.01 


+0.01 
-0.01 


0.02 


70, 100 


0.116 


± 0.008 


+0.018 
-0.003 


+0.005 
-0.005 


0.018 



Table 2: Values of the differential cross section as a function of Q 2 for the 
production of diffractive dijets. The range over which the cross section is averaged 
is given in the first column. Other details as in the caption of TableUl 



W bin 

(GeV) 


da/dW 
(pb/ GeV) 


"stat 

(pb/ GeV) 


^syst 

(pb/ GeV) 


<5es 

(pb/ GeV) 


Adiffr 
(pb/ GeV) 


100, 125 


0.26 


± 0.01 


+0.03 
-0.03 


+0.01 
-0.01 


0.01 


125, 150 


0.41 


± 0.02 


+0.04 
-0.03 


+0.02 
-0.03 


0.03 


150, 175 


0.67 


± 0.03 


+0.04 
-0.06 


+0.04 
-0.04 


0.03 


175, 200 


0.68 


± 0.02 


+0.03 
-0.04 


+0.05 
-0.04 


0.01 


200, 225 


0.77 


± 0.03 


+0.06 
-0.03 


+0.05 
-0.05 


0.05 


225, 250 


0.82 


± 0.03 


+0.03 
-0.06 


+0.05 
-0.05 


0.02 



Table 3: Values of the differential cross section as a function ofW. Other details 
as in the caption of Table [H 
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M x bin 
(GeV) 


da/dMx 
(pb/ GeV) 


"stat 

(pb/ GeV) 


Osyst 

(pb/ GeV) 


^ES 

(pb/ GeV) 


Adiffr 
(pb/ GeV) 


9, 14 


1.31 


± 0.07 


+0.02 
-0.08 


+0.05 
-0.06 


-0.03 


14, 20 


4.3 


± 0.1 


+0.2 
-0.2 


+0.2 
-0.2 


0.1 


20, 26 


4.5 


± 0.1 


+0.2 
-0.4 


+0.2 
-0.2 


0.0 


26, 32 


3.1 


± 0.1 


+0.2 
-0.3 


+0.3 
-0.2 


-0.1 


32, 42 


1.13 


± 0.05 


+0.08 
-0.06 


+0.12 
-0.09 


0.07 



Table 4: Values of the differential cross sections with respect to Mx- Other details 
as in the caption of Tabled 



(5 bin 


da/d/3 




O s tat 


"syst 


0~ES 


Adiffr 


(x 10" 2 ) 


(Pb) 




(pb) 


( P b) 


( P b) 


(Pb) 


0.32, 0.63 


1220 


± 


102 


+30 
-75 


+148 
-135 


-69 


0.63, 1.26 


2124 


± 


94 


+153 
-221 


+196 
-177 


-11 


1.26, 2.51 


1736 


± 


62 


+108 
-133 


+112 
-109 


46 


2.51, 5.01 


923 


± 


32 


+40 

-83 


+55 
-50 


3 


5.01, 10.00 


324 


± 


12 


+9 

-18 


+14 
-17 


3 


10.00, 19.95 


81.8 


± 


3.8 


+4.3 
-2.7 


+3.5 
-4.1 


4.1 


19.95, 39.81 


9.7 


± 


0.8 


+0.5 
-0.5 


+0.5 
-0.6 


0.4 



Table 5: Values of the differential cross sections with respect to (3. Other details 
as in the caption of Tabled 
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x°p S bin 
(x 1(T 2 ) 


xf s da/dxf s 
(pb) 


"stat 
( P b) 


"syst 
( P b) 


5es 
(Pb) 


Adiffr 
(Pb) 


0.25, 0.50 


24.3 


± 


1.8 


+0.8 
-1.5 


+1.0 

-l.i 


-0.5 


0.50, 0.79 


93 


± 


5 


+1 

-1 


+4 
-5 





0.79, 1.26 


195 


± 


7 


+3 
-9 


+9 
-10 


2 


1.26, 1.99 


306 


± 


10 


+10 
-25 


+17 
-17 


1 


1.99, 3.00 


409 


± 


13 


+33 
-33 


+35 
-30 


26 



Table 6: Values of the differential cross sections with respect to x^. Other details 
as in the caption of Tabled 



j bin 

(GeV) 


da/dE^ j 
(pb/ GeV) 


"stat 

(pb/ GeV) 


"syst 

(pb/ GeV) 


0~ES 

(pb/ GeV) 


Adiffr 
(pb/ GeV) 


4, 5.5 


51.7 


± 1.4 


+3.3 
-3.7 


+2.6 
-2.9 


2.9 


5.5, 7.5 


39.8 


± 1.1 


+2.6 
-2.8 


+2.2 
-2.0 


1.8 


7.5, 9.5 


9.7 


± 0.3 


+0.7 
-0.9 


+0.8 
-0.9 


0.2 


9.5, 11.5 


2.3 


± 0.1 


+0.1 
-0.1 


+0.3 
-0.2 


0.1 


11.5, 13.5 


0.65 


± 0.06 


+0.03 
-0.01 


+0.08 
-0.11 


0.03 


13.5, 16 


0.11 


± 0.02 


+0.02 
-0.02 


+0.01 
-0.03 


0.00 



Table 7: Values of the differential cross sections with respect to E^j- Other 
details as in the caption of Tabled 
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ti% bin 


da /dr]j 
(pb) 


^stat 

( P b) 


r 

^syst 

( P b) 


5pq 

(pb) 


/Xt^tddd 
' — l JJlr r ix 

(pb) 


-3.5, -3.0 


56.6 


± 1.9 


+2.5 
-3.8 


+7.8 
-7.6 


1.6 


-3.0, -2.5 


98.8 


± 2.9 


+3.6 
-6.2 


+7.2 
-7.1 


1.8 


-2.5, -2.0 


89.6 


± 2.6 


+5.7 
-6.0 


+5.1 
-4.9 


4.8 


-2.0,-1.5 


66.1 


± 2.1 


+4.1 
-4.2 


+3.7 
-4.1 


3.4 


-1.5, -1.0 


35.2 


± 1.3 


+3.3 
-2.6 


+2.7 
-2.0 


3.0 


-1.0, -0.5 


13.2 


± 0.7 


+1.4 
-1.3 


+1.3 
-1.3 


1.1 


-0.5, 0.0 


2.1 


± 0.2 


+0.4 
-0.5 


+0.4 
-0.3 


-0.2 



Table 8: Values of the differential cross sections with respect to rjj. Other details 
as in the caption of Tabled 



yobs 


da/dz°* s 
(Pb) 


"stat 
( P b) 


"syst 

(Pb) 


#ES 

(Pb) 


Adiffr 
(pb) 


0, 0.125 


24.5 


± 2.0 


+0.9 
-2.3 


+3.5 
-2.0 


-1.9 


0.125, 0.25 


134.6 


± 5.4 


+8.8 
-13.8 


+12.3 
-10.8 


0.7 


0.25, 0.375 


155.1 


± 5.7 


+9.6 
-12.5 


+10.1 
-9.8 


5.2 


0.375, 0.5 


133.7 


± 5.1 


+8.3 
-10.2 


+6.1 
-8.1 


5.7 


0.5, 0.625 


100.6 


± 4.2 


+5.8 
-7.8 


+5.2 
-5.4 


2.4 


0.625, 0.75 


80.4 


± 3.6 


+1.3 
-2.8 


+3.8 
-3.8 


0.5 


0.75, 0.875 


55.5 


± 2.8 


+1.7 
-3.1 


+2.8 
-2.8 


-1.5 


0.875, 1 


31.5 


± 2.1 


+3.1 
-4.0 


+2.2 
-1.5 


-1.3 



Table 9: Values of the differential cross sections with respect to . Other details 
as in the caption of Tabled 
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x° bs bin 


da/dxf s 
(Pb) 


Ogtat 

(Pb) 


"syst 

(Pb) 


^ES 
(pb) 


Adiffr 
(Pb) 


0, 0.25 


5.3 


± 0.6 


+0.8 
-0.5 


+0.6 
-0.2 


0.7 


0.25, 0.5 


25.0 


± 1.3 


+2.8 
-1.5 


+1.5 
-2.2 


2.6 


0.5, 0.75 


87.4 


± 3.1 


+4.5 
-5.7 


+7.8 
-8.2 


3.5 


0.75, 1 


240.5 


± 6.7 


+ 11.0 
-17.2 


+12.4 
-11.7 


5.9 



Table 10: Values of the differential cross sections with respect to x° . Other 
details as in the caption of Tabled 
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zf s bin 


da/dz^dE^ 
(pb/ GeV) 


$stat 

(pb/ GeV) 


"syst 

(pb/ GeV) 


<$ES 

(pb/ GeV) 


Adiffr 
(pb/ GeV) 


5.0 < E^ dl < 6.5GeV(< E* Tjl >= 5.7 GeV) 


0, 0.25 


31.9 


± 1.5 


+2.1 
-4.2 


+1.9 
-1.4 


-1.2 


0.25, 0.375 


53.1 


± 2.6 


+4.2 
-4.6 


+2.3 
-2.3 


3.0 


0.375, 0.5 


46.7 


± 2.4 


+3.0 
-3.0 


+1.9 
-2.4 


2.5 


0.5, 0.625 


35.3 


± 2.1 


+1.1 
-1.7 


+1.6 
-1.8 


1.0 


0.625, 0.75 


29.3 


± 1.9 


+0.3 
-1.9 


+1.3 
-1.2 


-0.8 


0.75, 0.875 


18.4 


± 1.4 


+1.2 
-1.9 


+0.8 
-0.8 


-1.0 


0.875, 1 


11.4 


± 1.2 


+0.3 
-1.0 


+0.5 
-0.5 


-0.5 


6.5 < < 8.0GeV(< E^ x >= 7.2 GeV) 


0, 0.25 


13.2 


± 0.5 


+1.4 

-1.5 


+0.9 
-0.6 


0.17 


0, 0.25 


13.2 


± 0.8 


+1.4 
-1.5 


+0.9 
-0.6 


0.2 


0.25, 0.375 


25.9 


± 1.5 


+1.2 
-2.2 


+1.9 
-2.0 


-0.9 


0.375, 0.5 


21.9 


± 1.3 


+1.9 
-1.8 


+1.6 
-0.9 


1.9 


0.5, 0.625 


18.3 


± 1.2 


+0.8 
-1.0 


+1.2 
-1.0 


0.7 


0.625, 0.75 


14.8 


± 1.1 


+0.8 
-0.9 


+0.6 
-0.6 


0.2 


0.75, 0.875 


12.4 


± 1.0 


+0.8 
-0.9 


+0.6 
-0.8 


0.0 


0.875, 1 


5.6 


± 0.7 


+0.2 
-0.2 


+0.5 
-0.2 


-0.1 


8.0 < Etpjx < 16.0GeV(< E^ jX >= 9.7 GeV) 


0.25, 0.375 


2.4 


± 0.2 


+0.4 
-0.3 


+0.2 
-0.2 


0.2 


0.375, 0.5 


2.4 


± 0.2 


+0.2 
-0.2 


+0.2 
-0.3 


-0.1 


0.5, 0.625 


1.9 


± 0.1 


+0.1 
-0.2 


+0.2 
-0.2 


0.0 


0.625, 0.75 


1.7 


± 0.1 


+0.1 
-0.0 


+0.2 
-0.2 


0.1 


0.75, 0.875 


1.4 


± 0.1 


+0.0 
-0.1 


+0.1 
-0.1 


0.0 


0.875, 1 


0.80 


± 0.09 


+0.01 
-0.10 


+0.10 
-0.08 


-0.05 



Table 11: Values of the double differential cross sections with respect to z^ in 
bins of -Ej-ji- Other details as in the caption of Tabled 
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4 bs bin 


da/dzf s dQ 2 
(pb/ GeV 2 ) 


"stat 

(pb/ GeV 2 ) 


"syst 

(pb/ GeV 2 ) 


o~es 

(pb/ GeV 2 ) 


Adiffr 
(pb/ GeV 2 ) 


5 < Q 2 < 12 GeV 2 (< Q 2 >= 8.1 GeV 2 ) 


0, 0.25 


5.1 


± 0.3 


+0.4 
-0.5 


+0.1 
-0.4 


0.0 


0.25, 0.375 


8.7 


± 0.4 


+0.6 
-0.7 


+0.7 
-0.7 


-0.1 


0.375, 0.5 


7.2 


± 0.4 


+0.4 
-0.5 


+0.5 
-0.5 


0.3 


0.5, 0.625 


5.2 


± 0.3 


+0.1 
-0.2 


+0.4 
-0.4 


0.1 


0.625, 0.75 


4.3 


± 0.3 


+0.1 
-0.2 


+0.3 
-0.2 


0.1 


0.75, 0.875 


2.9 


± 0.2 


+0.1 
-0.1 


+0.2 
-0.2 


-0.1 


0.875, 1 


1.5 


± 0.2 


+0.1 
-0.2 


+0.1 
-0.1 


-0.1 


12 < Q 2 < 25 GeV 2 (< Q 2 >= 17.2 GeV 2 ) 


0, 0.25 


1.43 


± 0.09 


+0.14 
-0.16 


+0.09 
-0.07 


0.08 


0.25, 0.375 


3.0 


± 0.2 


+0.4 
-0.2 


+0.1 
-0.2 


0.4 


0.375, 0.5 


2.3 


± 0.1 


+0.1 
-0.2 


+0.1 
-0.1 


0.1 


0.5, 0.625 


2.0 


± 0.1 


+0.1 
-0.2 


+0.1 
-0.1 


0.0 


0.625, 0.75 


1.6 


± 0.1 


+0.1 
-0.1 


+0.1 
-0.1 


0.0 


0.75, 0.875 


1.2 


± 0.1 


+0.0 
-0.1 


+0.1 
-0.1 


-0.1 


0.875, 1 


0.61 


± 0.07 


+0.01 
-0.03 


+0.04 
-0.03 


-0.01 



Table 12: Values of the double differential cross sections with respect to z°p in 
bins of -Erji- Other details as in the caption of Tabled 
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z° bs bin 


da/dzf s dQ 2 
(pb/ GeV 2 ) 


"stat 

(pb/ GeV 2 ) 


"syst 

(pb/ GeV 2 ) 


5es 

(pb/ GeV 2 ) 


Adiffr 
(pb/GeV 2 ) 


25 < Q 2 < 50 GeV 2 (< Q 2 >= 35.2 GeV 2 ) 


0, 0.25 


0.51 


± 0.04 


+0.03 
-0.08 


+0.03 
-0.02 


-0.06 


0.25, 0.375 


1.03 


± 0.07 


+0.07 
-0.13 


+0.06 
-0.07 


-0.03 


0.375, 0.5 


1.00 


± 0.07 


+0.06 
-0.06 


+0.04 
-0.06 


0.05 


0.5, 0.625 


0.77 


± 0.06 


+0.06 
-0.02 


+0.04 
-0.04 


0.06 


0.625, 0.75 


0.60 


± 0.05 


+0.03 
-0.04 


+0.04 
-0.03 


-0.01 


0.75, 0.875 


0.44 


± 0.04 


+0.01 
-0.02 


+0.03 
-0.04 


0.00 


0.875, 1 


0.24 


± 0.03 


+0.01 
-0.01 


+0.02 
-0.01 


0.00 


50 < Q 2 < 100 GeV 2 (< Q 2 >= 69.5 GeV 2 ) 


0, 0.25 


0.10 


± 0.01 


+0.00 
-0.01 


+0.01 
-0.01 


0.00 


0.25, 0.375 


0.25 


± 0.02 


+0.05 
-0.01 


+0.02 
-0.01 


0.05 


0.375, 0.5 


0.28 


± 0.03 


+0.02 
-0.01 


+0.01 
-0.02 


0.02 


0.5, 0.625 


0.20 


± 0.02 


+0.02 
-0.01 


+0.01 
-0.01 


0.02 


0.625, 0.75 


0.16 


± 0.02 


+0.01 
-0.00 


+0.01 
-0.01 


0.01 


0.75, 0.875 


0.13 


± 0.02 


+0.00 
-0.00 


+0.01 
-0.01 


0.00 


0.875, 1 


0.11 


± 0.02 


+0.01 
-0.02 


+0.01 
-0.01 


-0.01 



Table 13: Values of the double differential cross sections with respect to in 
bins of E^ii- Other details as in the caption of Tabled 
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Figure 2: The measured 7?max distribution (dots) (a) before diffractive selection, 
(b) after the Eppc cut and (c) after adding the x°^ s cut. Also shown are area- 
normalised MC expectations obtained by fitting the relative amount o/Rapgap and 
Djangoh to give the best description of the data before any diffractive selection. 
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Figure 3: Measured differential cross section (dots) as a function of (a) Q 2 , (b) 
W, (c) Mx, (d) (3 and (e) x^. The inner error bars represent the statistical uncer- 
tainty and the outer error bars represent the statistical and uncorrelated systematic 
uncertainties added in quadrature. The shaded band represents the correlated un- 
certainty. For comparison the area-normalised predictions of the Rapgap (solid 
lines) and the Satrap (dashed lines) MC models are also shown. 
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Figure 4: Measured differential cross section as a function of (a) J; (b) rjj, (c) 
z^ s and (d) x° bs . The dashed- dotted line represents the area-normalised Rapgap 
with only the direct photon contribution. Other details as in the caption of Fig. 
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Figure 5: Measured differential cross section as a function of (a) Q 2 , (b) W , 
(c) Mx, (d) (3 and (e) x°p compared to the NLO predictions obtained using the 
available dPDFs, as indicated in the figure. The hatched area indicates the theoret- 
ical uncertainty of the predictions estimated using the ZEUS LPS+charm dPDFs. 
Other details as in the caption of Fig. [3J 
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Figure 6: Measured differential cross section as a function of (a) E^ 3 , (b) rjj, 
(c) z°p S and (d) x° bs compared to the NLO prediction obtained using the available 
dPDFs. Other details as in the caption of Fig. 
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Figure 7: Ratio, R, of the data to the NLO prediction using the ZEUS LPS+charm 
dPDFs (dots) as function of (a) Q 2 , (b) W , (c) M x , (d) (3 and (e) x°p S . Also shown 
is the ratio of NLO calculations with other dPDFs to ZEUS LPS+charm. Other 
details as in the caption of Fig. 
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Figure 8: Ratio, R, of the data to the NLO prediction using the ZEUS LPS+charm 
dPDFs (dots) as function of (a) E^ 3 , (b) r]j, (c) z° bs and (d) x° hs . Other details 
as in the caption of Fig. \7\ 
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Figure 9: Measured differential cross section as a function of z^ s in different 
regions of • 1 (dots). Other details as in the caption of Fig. 
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Figure 10: Measured differential cross section as a function of z^ s in different 
regions of Q 2 (dots). Other details as in the caption of Fig. [5j 
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